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Abstract. We use lattice QCD to calculate the B-mixing hadronic matrix elements for a basis of 
effective four-quark operators that spans the space of all possible contributions in, and beyond, 
Q-i the Standard Model. We present results for the SU(3)-breaking ratio | and discuss our ongoing 

calculation of the mixing matrix elements, including the first calculation of the beyond the Standard 
Model matrix elements from unquenched lattice QCD. 
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a MOTIVATION 
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Through a combination of the GIM mechanism, Cabibbo suppression, and loop suppres- 
sion, the Standard Model (SM) contribution to 5-mixing is small and new physics (NP) 
effects may be discernible [1]. In fact, there are experimental hints this may be the case. 
In unitarity triangle analyses [2, 3] a persistent 2 — 3(7 inconsistency is suggestive of NP 
and points to 5-mixing as a possible source. D0's same-sign dimuon charge asymme- 
try [4] and global analyses by UTfit [5] and CKMfitter [6] reveal 2 — 4a discrepancies 
in SM 5-mixing. Experimental measurements of the 5-mixing oscillation frequency [7] 
have sub-percent precision but cannot be fully leveraged in the search for NP as theory 
errors, dominated by hadronic uncertainty, are an order of magnitude larger [8]. 
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CALCULATION 

To lowest order in the SM, 5-mixing is described by box diagrams [cf. Fig. 1 (left)]. Un- 
der the operator product expansion (OPE) flavor-changing short-distance interactions, of 
0(100 GeV) in the SM and higher in NP scenarios, and long-distance hadronic physics, 
of 0(500 MeV), factorize. At the energies relevant to 5-mixing, of 0(A/g ~ 5 GeV), 
flavor-changing physics is described by the local, effective, four-quark interaction of 
Fig. 1 (right). A commonly used basis of mixing operators is 

ff x = (b^Lq") (b^Lq P ) ^ = (b a Lq a ) (b^RqP) 

ff 2 = {b a Lq a ) (bVLqV) 5 = (b a Lq^) (b P Rq a ) 

3 = {b a Lqt)(pLq a ) (1) 



where L/R are left/right projection operators and a,/3 are color indices. A calculation 
of the matrix elements of these operators is sufficient to parameterize the hadronic 




FIGURE 1. (Left) A SM contribution to B-mixing and (right) an effective four-quark interaction. 



contributions to 5-mixing in and beyond the SM. 

Under the OPE, the expression for the oscillation frequency factorizes: 

AM, = £Q(Ai) (B° q \^)\B° q ), (2) 

where the short-distance C, are model dependent and the long-distance hadronic mixing 
matrix elements (B q \&i\B q ) must be calculated nonperturbatively using lattice QCD. 

The phenomenologically-relevant bag parameters B$ and SU(3)-breaking ratio £ are 
related to the matrix elements by 

(B q \^\B° q ) = afl/l and % = f Bs ^/f Bd ^, (3) 

where fs is the 5-meson decay constant and q are numerical factors. 1 

The Lattice Calculation. Working in the 5-meson rest-frame, we generate correla- 
tion function data via Monte Carlo evaluation of the path integral representations of the 
vacuum expectation values 

C 2 ?\t) = (B q (t)B q (0f) and C^(t h t 2 ) = (B Q q {t 2 )^{Q)B Q q { h )) : (4) 

with gauge field integration performed with the MILC ensembles. We simulate with 
staggered light and Fermilab heavy quarks (for details of gluon and quark discretizations 
see [9] and references therein). Using Bayesian fitting techniques [10], these data are fit 
to the ansatze [11], 

Cf pt (t u t 2 ) = £ (B° q \0i\B q ) ZnZm e -Mnh-M m t 2 ^ (5) 

n,m=0 2y/M n M m 

where Z n is the amplitude and M n the mass of the 5-meson n th excited state and T 
the temporal extent of the lattice. From these fits we extract the matrix elements over 



An alternate definition of the bag parameter includes a factor of M\ / (mj + m q ) 2 with c,- for i ^ 1. 



a range of valence-quark masses, light sea-quark masses, and lattice spacings. Lattice 
matrix elements are matched to the continuum at one loop in tadpole-improved lattice 
perturbation theory [12]. The continuum matrix elements corresponding to physical light 
(for B Q d ) or strange (for B®) valence quark, and at physical light sea-quark mass, are 
then obtained by extrapolation/interpolation with the aid of rooted, staggered, chiral 
perturbation theory [13, 14]. 

Status and Outlook. We recently completed a calculation of t, = 1.268(63) using 
data at lattice spacings down to ~ 0.09 fm and valence quarks as light as 0. 1 m s [15]. 

An ongoing calculation includes an update of ^ and the calculation of matrix elements 
and bag parameters for all five operators in Eq. (1). In addition, it includes several 
improvements: a three-fold increase in statistics; data at lattice spacings as small as 
~ 0.045 fm and valence quarks as light as 0.05 m s , to reduce the effect of the continuum- 
chiral extrapolation; and a more thorough treatment of chiral perturbation theory [14]. 

Preliminary results for the matrix elements with the new, extended data set can be 
found in Table 4 of Ref. [16]. These results are based on data at lattice spacings down 
to f^O.06 fm and valence quarks as light as 0.1 m s . Initial studies show that our new 
treatment of chiral perturbation theory has a negligible effect on the values of the matrix 
elements in this analysis. 
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